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ABSTRACT 

We present refined values for the physical parameters of transiting exoplanets, based on a self- 
consistent and uniform analysis of transit light curves and the observable properties of the host stars. 
Previously it has been difficult to interpret the ensemble properties of transiting exoplanets, because of 
the widely different methodologies that have been applied in individual cases. Furthermore, previous 
studies often ignored an important constraint on the mean stellar density that can be derived directly 
from the light curve. The main contributions of this work are i) a critical compilation and error 
assessment of all reported values for the effective temperature and metallicity of the host stars; ii) 
the application of a consistent methodology and treatment of errors in modeling the transit light 
curves; and Hi) more accurate estimates of the stellar mass and radius based on stellar evolution 
models, incorporating the photometric constraint on the stellar density. We use our results to revisit 
some previously proposed patterns and correlations within the ensemble. We confirm the mass-period 
correlation, and we find evidence for a new pattern within the scatter about this correlation: planets 
around metal-poor stars are more massive than those around metal-rich stars at a given orbital period. 
Likewise, we confirm the proposed dichotomy of planets according to their Safronov number, and we 
find evidence that the systems with small Safronov numbers are more metal-rich on average. Finally, 
we confirm the trend that led to the suggestion that higher-metallicity stars harbor planets with a 
greater heavy-element content. 

Subject headings: methods: data analysis — planetary systems — stars: abundances — stars: fun- 
damental parameters — techniques: spectroscopic 



1. INTRODUCTION 

The transiting exoplanets are only a small subset 
of all the known planets orbiting other stars, but 
they hold tremendous promise for deepening our un- 
derstanding of planetary formation, structure, and evo- 
lution. Observations of transits and occultations 3 
(along with the spectroscopic orbit of the host star) 
not only allow one to measure the mass and ra- 
dius of the planet, but also provide op portunities to 
meas u re the stella r spin-o rbit alignment (jQueloz et alJ 
120001 I Winn et alJ l2007al). the planet a ry brightness 
temp erature (|Charbonneau et all 120051 : iDeming et al.l 
I2005D. the planetary day-night temperature difference 
(jKnutson et all2007j ). and eve n absorption lines of plane- 
tary atmospheric constituents (ICharbonneau et al. 2002; 
IVidal-Madiar et al.|[200llTinetti et al II2007D . These and 
other observations have been accompanied by theoretical 
progress in modeling the physical processes in the plan- 
etary interiors and atmospheres, as well as the planets' 
interactions with their parent stars. This rapid progress 
has been stimulated in no small measure by the remark- 
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3 The word transit is sometimes assumed in the field of exoplanet 
research to be synonymous with eclipse. In reality, it has a more 
restricted meaning and has long been used in the eclipsing binary 
field to describe an eclipse of the larger object by the smaller one. 
The term occultation is used to refer to the passage of the smaller 
object (in this case, the planet) behind the larger one (the star) 
(see, e.g., [Popper 197rJ). To avoid confusion, we advocate that 
occultation or secondary eclipse are preferable to neologisms such 
as "secondary transit" or "anti-transit". 



able diversity of planet characteristics that has been 
found among the members of the transiting ensemble. 

The accuracy and precision with which the properties 
of the planet can be derived from transit data depend 
strongly on whatever measurements and assumptions are 
made regarding the host star. For example, for a given 
value of the transit depth, the inferred planetary radius 
scales in proportion to the assumed stellar radius; and for 
a given spectroscopic orbit of the host star, the inferred 
planetary mass scales as the two-thirds power of the stel- 
lar mass. In the literature on transiting planets, a wide 
variety of methods have been used to estimate the radius 
and mass of the parent star, ranging from simply looking 
them up in a table of average stellar properties as a func- 
tion of spectral type, all the way to fitting detailed stel- 
lar evolutionary models constrained by the luminosity, 
effective temperature, and other observations that may 
be available for the star. As a result, the ensemble of 
planet properties at our disposal is inhomogeneous, and 
in many cases the uncertainty of those determinations is 
dominated by systematic errors in the stellar parameters 
that are treated differently by different investigators. 

This situation is unfortunate because it hinders our 
ability to gauge the reliability of any patterns that are 
discerned among the ensemble properties of transiting 
exoplanets. With 23 systems that have been reported 
in the literature, this subfield should be poised for the 
transition from a handful of results to a large and diverse 
enough sample for meaningful general conclusions to be 
drawn, but the heterogeneity of reported results is clearly 
an obstacle. 

It may be surprising that we are limited in many cases 
by our knowledge of the properties of the parent stars; 
one would think that stellar physics is a "solved prob- 
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lem" in comparison to exoplanetary physics. However 
it must be remembered that the host stars are usually 
isolated (and therefore no dynamical mass measurement 
is possible), and that many of the host stars are dis- 
tant enough that they do not even have measured par- 
allaxes. Of the 23 cases in the literature , five have Hip- 
parcos parallaxes (|Perrvman et al.ll"l997( ). For those few 
it is fairly straightforward to estimate the stellar prop- 
erties, but for the other systems, more indirect methods 
have been used. These indirect methods often rely on 
the value of the stellar surface gravity that is derived 
by measuring the depths and shapes of gravity-sensitive 
absorption lines in the stellar spectrum. This is a no- 
toriously difficult measurement and the result is often 
strongly correlated with othe r parameters that affe ct the 
spectrum. Recently, however. iSozzetti et al.l (|2007l ). Hol- 
man et al. (2007), and others demonstrated that it is 
possible to do better by using the information about the 
mean stellar density that is encoded in the transit light 
curve. This information was typically overlooked prior 
to these studies. 

The study presented in this paper was motivated by the 
desire for a homogeneous analysis, and by the desire to 
take advantage of the photometric estimates of the stellar 
mean density. We have revisited the determination of the 
stellar parameters for all of the transiting planets that 
have been reported in the literature. We have taken the 
opportunity to merge all existing measurements of the 
atmospheric parameters (mainly the effective tempera- 
ture and metallicity) with the goal of presenting the best 
possible values. We have chosen a uniform method for 
analyzing photometric data and have re-analyzed exist- 
ing light curves where necessary to provide homogeneity. 
Our hope was that by applying these procedures across 
the board, we and other investigators could view the en- 
semble properties with greater clarity and uncover any 
interesting clues the transiting planets might provide us 
about the origin, structure, and evolution of exoplanets. 

This paper is organized as follows. §[2] describes the 
procedures by which we estimated the stellar proper- 
ties, using the available spectroscopic and photometric 
datasets, and a particular set of theoretical stellar evolu- 
tion models. As a check on the models, §[3] compares the 
results of a subset of our calculations with those derived 
from a different set of evolutionary models. §2] investi- 
gates alternate ways of estimating the stellar properties. 
§[5]deals with GJ 436, which needs special treatment be- 
cause the host star has such a lower mass than the other 
host stars. §[6] presents the final results for the planetary 
parameters. §[7] uses the new results to check on some 
of the previously proposed correlations among the prop- 
erties of the transiting ensemble, and §|5] provides final 
remarks. The Appendix lists the data sets and other 
issues that are particular to each system. 

2. DETERMINING THE STELLAR PROPERTIES 

For a transiting planet, the basic data are the spectro- 
scopic orbit of the star (radial- velocity curve), and the 
photometric observations of transits (light curve) . With 
such data, the planetary mass and radius cannot be de- 
termined independently of the stellar properties. The 
radial-velocity curve can be used to determine 

M p smi = 4.9f9xl0" 3 P 1/3 (l-e 2 ) 1/2 ^ [(Af* + M p )/M Q ] 2/3 

(1) 



(in units of the mass of Jupiter) , where i is the inclination 
angle of the orbit, P is the orbital period in days, e is 
the eccentricity, and K+ the velocity semi-amplitude of 
the star in ms . Even when sini is known precisely, 
the value of the planetary mass M p that is derived from 

the data will scale as M* ' , where M* is the mass of the 
star. Meanwhile, the light curve does not immediately 
yield R p , the planetary radius; rather, the ratio of the 
radii R p / ' R* is pinned down through the depth of the 
transit. 

The stellar mass and radius are usually inferred in- 
directly, from an analysis of high-resolution spectra of 
the star with the aid of model atmospheres, followed by 
a comparison of the atmospheric parameters with stellar 
evolution models. The latter step is performed somewhat 
differently by different authors depending on the obser- 
vational constraints available. For this work we have 
compiled and critically reviewed all of the available in- 
formation regarding the atmospheric properties of the 
host stars. This effort is described in § 12.11 

The approach adopted in this paper makes use of in- 
formation from the light curves of transiting planets to 
constrain the mean density of the star, p+, following 
ISozzetti et al.l f|2007T ) . As described therein, the quantity 
a/Ri, (the planet-star separation a in units of the stel- 
lar radius) is directly related to the mean stellar density 
(hereafter, simply the density), and can be derived from 
the photometry without much knowledge about the star 
(see below). The only dependence o/i?* has on the stel- 
lar properties is through the limb-darkening coefficients, 
which is typically a second-order effect (see § I2.2[) , 

Many of the published light curve analyses do not re- 
port the value of a/R+ explicitly, and it is difficult or 
impossible to reconstruct its value accurately from the 
published information, ft is especially difficult to obtain 
a good measure of the true uncertainty in a/R+ from 
the published information. For this reason, and for the 
sake of homogeneity, we have re-analyzed many of the 
high-quality photometric time-series available to us for 
all transiting planets. We describe this effort in § 12.21 

2.1. Atmospheric parameters 

High-resolution spectroscopic studies have been under- 
taken for almost all of the parent stars of the known 
transiting planets. The basic products of these studies 
are measurements of the effective temperature T c g, iron 
abundance [Fe/H], and surface gravity log g. The quoted 
precision of these determinations varies widely, depend- 
ing on the signal-to-noise ratio and resolution of the ob- 
servations, the modeling techniques employed, and the 
attitude taken toward systematic errors. 

For this work we have relied on published determina- 
tions, rather than any new spectroscopic data. In many 
cases, a given transiting system has been described by 
more than one analysis, and they do not necessarily 
agree. In those cases, rather than arbitrarily choosing 
the most recent study, or the one with the smallest uncer- 
tainties, we critically examined all of the available studies 
and combined them. Our choices are documented in the 
Appendix. We were guided by our own experience and 
increased the error estimates whenever they seemed op- 
timistic. For example, many of the automated spectro- 
scopic analysis tools in common use today return formal 
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uncertainties for the effective temperatures that are only 
a few tens of degrees Kelvin; for one of the transiting 
systems the quoted error was only 13 K. There is ample 
literature on the subject of the absolute effective temper- 
ature scale and the systematic errors inherent in placing 
any si ngle system on such a scale . A recent investiga- 
tion bv lRarmrez fc M elcndcz ( 2003) has shown that there 
are still differences of order 100 K between temperatures 
derived from the spectroscopic condition of excitation 
equilibrium and from the Infrared Flux Method (IRFM, 
IBlackwell fc Shallid [l977t IBlackwell eTail [l980h . Other 
studies have indicated discre pancies of 50-100 K between 
different temperature sc ales (| Ramirez fc Melendezl l2005: 
iCasagrande et al.l [2006) , and discussed at length possi- 
ble sources for these errors. In light of these findings, 
we considered it prudent to adopt temperature errors no 
smaller than 50 K for our study, and then only when 
there are several independent and consistent determina- 
tions or other evidence supporting that level of accu- 
racy. Similar concerns hold for the spectroscopic sur- 
face gravities, although we do not actually make use of 
them in this work, as described below. For the spec- 
troscopic metallicities, we adopted a minimum uncer- 
tainty of 0.05 dex (for the best cases with multiple in- 
dependent measurements), even though smaller errors 
have occasionally been reported for individual analy- 
ses. This is mainly because of the stro ng correlations 
prese nt among [Fe/H], T c g, and logg (e.g.. lBuzzoni et ahl 
2001), as well as some evidence for sy stematic differences 
between different groups (see, e.g.. iSantos et al.l 12004 
iFischer fc Valentill2005l I Gonzalez fc Lawsll2007l l In the 
Appendix we provide a complete listing of the sources we 
have drawn from in each case. The values of T e s, [Fe/H], 
and log g finally adopted for all systems are listed in Ta- 
ble [T] Some degree of non-uniformity in these quantities 
is unavoidable due to the variety of procedures used by 
different authors, but we believe they represent the best 
available set for the parent stars based on current knowl- 
edge. 

2.2. Light curve fits 

For each system, we examined the highest-quality tran- 
sit photometry available to us in order to determine the 
key parameters R p /R+, a/R+, and i. Our methodol- 
ogy is described below, and the details of the data sets 
that were used in each case are given in the Appendix. 
In some cases, the published determinations of R p /R+, 
a/R*, and i matched our own methodology very closely, 
and we simply adopted the values from the literature; 
these cases are also specified in the Appendix. 

In the absence of limb darkening, the four primary ob- 
servables in a transit light curve are the midtransit time, 
the depth, the total duration, and the partial-phase dura- 
tion (ingress or egress). A sequence of measured midtran- 
sit times usually leads to a very precise determination of 
the orbital period P. The depth is equal to the planet-to- 
star radius ratio, (R p / R±) 2 . At fixed P, the parameters i 
and a/R* can be written in terms of the total and partial 
dur ations through an application o f Kepler's Law (see, 
e.g. JSeager fc Mallen-Ornelasl[2003l ). With limb darken- 
ing, however, there is no longer a well-defined depth or 
partial duration; an accurate light curve model must in- 
clude a realistic intensity distribution across the stellar 
disk, which will depend on the stellar temperature, sur- 



face gravity, and metallicity. Thus, to some degree, the 
determination of a/i?*, and i must be accompa- 

nied by some assumptions about the stellar properties. 

Our procedure was as follows. We modeled each sys- 
tem using a two-body Keplerian orbit. The star has mass 
M+ and radius i?*, and the planet has mass M p and ra- 
dius R p . The orbit has period P, eccentricity e, argument 
of pericenter to, and inclination i. For our purpose the 
uncertainty in P was completely negligible; we fixed P at 
the most precisely determined value in the literature. In 
almost all cases, the radial-velocity data are consistent 
with a circular orbit, and we assumed e = exactly. For 
HAT-P-2 (HD 147506) we fixed the values of e and u> at 
those that have been derived from radial-velocity data 
(and in the end we verified that changing these param- 
eters by la does not significantly affect any of our final 
results). The initial condition is specified by a particu- 
lar midtransit time T c . When the sky projections of the 
star and planet do not overlap, the model flux is unity. 
When they do overlap we use the analytic formulas of 
iMandel fc Ago] (|2002f) to compute the integral of the in- 
tensity over the unobscured portion of the stellar disk, 
assuming a quadratic limb-darkening law. 

To arrive at a self-consistent solution including limb 
darkening, we began with initial values for M* and M p 
from the literature. We also chose values for the stel- 
lar T e ff, log<?, and metallicity, as described in the previ- 
ous section. Then we adopted limb-darkening coefficients 
based on those s tellar param eters, by interpolating the 
tables of IClaretl (|2000L \2004 ) for the appropriate band- 
pass. At this point we estimated all of the remaining 
parameters (using the procedure described in the next 
paragraph) and used the result for a/i?* to update the 
determination of the stellar properties. Then, a new pho- 
tometric parameter estimation was performed, with re- 
vised values of M+, M p and the limb-darkening coeffi- 
cients, and so forth. This procedure converged after two 
or three iterations, in the sense that further iterations 
changed none of the parameters by more than about a 
tenth of the statistical error. 

The parameter estimations were carried out with a 
Mar kov Chain Mon t e Ca rlo algorithm (MCMC; see, 
e.g., iTegm ark et al. (2004) for applications to cosmo- 
logical data, Fordl (120051) for radial - veloc ity data, and 



iHolman et al.l (|2006fl or lBurke et~aTI (|2007l ) for a similar 
approach to transit fitting) . It is based on the goodness- 
of-fit statistic 
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/j(obs) - ^ (calc) 



(2) 



where fj (obs) is the flux observed at time j and o~j con- 
trols the weights of the data points, and /j(calc) is the 
flux calculated with our model. In the MCMC algorithm, 
a stochastic process is used to create a sequence of points 
in parameter space whose density approximates the joint 
a posteriori probability density for all parameters. One 
begins with an initial point and iterates a jump func- 
tion, which in our case was the addition of a Gaussian 
random number ("perturbation") to a randomly chosen 
parameter. If the new point has a lower x 2 , the jump is 
executed; if not, the jump is executed with probability 
exp(— A% 2 /2). We adjusted the sizes of the perturba- 
tions until approximately ~25% of jumps are executed 
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for each parameter. 

For the data weights <jj , we used the observed standard 
deviation of the out-of-transit data (ci), multiplied by a 
factor (3 > 1 to account at least approximately for time- 
correlated errors ( "red" noise) which are often significant 
for ground-based data. We chose (3 as follows. The key 
timescale is the partial-phase duration, because the lim- 
iting error in the determination of a/R* and i is generally 
the fractional error in the partial-phase duration. We av- 
eraged the out-of-transit data over this timescale, with 
each time bin consisting of N points depending on the 
cadence of observations, and then calculated the stan- 
dard deviation of the binned data, o~n- Finally we set 
P = un * VN jo\. With white noise only, we would ob- 
serve P — 1, but in practice P > 1 because the number 
of effectively independent data points is smaller than the 
actual number of data points. For this paper we deliber- 
ately chose to analyze only those data for which P < 2. 

For each parameter, we took the mode of the MCMC 
distribution after marginalizing over all other parameters 
to be the "best value." We defined the 68% confidence 
limits pio and phi as the values for which the integral of 
the distribution between p\ a and phi is 0.68, and the in- 
tegrals from the minimum value to p\ and from phi to 
the maximum value were both 0.16 = (1.00 — 0.68)/2. In 
some cases the mode was all the way at one end of the 
probability distribution; in particular there were several 
cases in which i = 90° was the mode of the inclination 
distribution. In those cases we report phi as 90° and p\ Q 
as the value for which the integral from zero to p\ Q was 
(1.00 - 0.68) = 0.32. The final results are given in Ta- 
ble [2] including the stellar density p+ computed directly 
from a/R-t and the period. 

2.3. Stellar masses, radii, luminosities, surface 
gravities, and ages 

The fundamental parameters for the host stars are de- 
rived here using stellar evolution models. We rely on the 
spectroscopically determined T e s and [Fe/H], and we re- 
quire also an indicator of luminosity (L*) or some other 
measure of evolution. Since most of these stars lack a 
parallax measurement, the spectroscopic surface gravity 
has often been used in the past as a luminosity indi- 
cator. The effect of log g on the spectral lines is rela- 
tively subtle, and strong correlations between logg and 
both temperature and metallicity make these deter mina- 
tions challenging. Sea ger fc Mallen-Ornelad (|2003h have 
pointed out that an important property intrinsic to the 
star, the density, is encoded in the transit light curves 
mainly through its dependence on the tr ansit duration. 
As sh own there, and more explicitly by Sozze tti et alJ 
(2007), the density is directly related to a/R*, one of the 
parameters often solved for in modeling the photometry. 
This quantity can typically be determined more precisely 
than log g and it is highly sensitive to the degree of evo- 
lution of the star, i.e., to its size. Thus, it serves as a 
better proxy for luminosity in many cases. We illustrate 
this below. 

To determine the stellar mass and radius, and other 
relevant properties of the host stars, we follow the proce- 
dure described by ISozzetti et al.l (|2007l ) with minor im- 
provements, and compare model isochroncs directly with 
the measured values of T c g, [Fe/H] , and a/R±. The latter 



can be calculated from the models as 

G_\ 1/3 P 2 / 3 
in 2 j R+ 



(3) 



where G is the Newtonian gravitational constant, and 
the period P is well known from the photometry. The 
planet mass M p is not known a priori, but its influence is 
very small compared to the stellar mass Af+, and even a 
rough value is usually sufficient for this application. Once 
the stellar mass is known, the process can be repeated if 
necessary with an improved value of M p . 

The stellar evolution model s we use are t hose from 
the Yonsei-Yale (Y 2 ) series bv lYi et all (pOOlh ( see also 
iDemarque et al.ll2004h . which are conveniently provided 
with tools for interpolating isochrones in both age and 
metallicity. 4 We explore the full range of metallicities 
allowed by the observational errors in [Fe/H], sampling 
~20 equally spaced values for each system. For each 
metallicity we consider a range of ages from 0.1 to 14 
Gyr, in steps of 0.1 Gyr. These isochrones are inter- 
polated to a fine grid in mass, and compared point by 
point with the measured values of T c g and a/i?*. All lo- 
cations ( "matches" ) on the isochrone that are consistent 
with these quantities within the observational errors are 
recorded. We also record the corresponding likelihood 
given by exp(— x 2 /2), where 



X 



/ A [Fe/H] \ | /AT eg 



Va/R, 



and the A quantities represent the difference between 
the observed and model values at each point. Observa- 
tional errors are assumed to be Gaussian, and the asym- 
metric error bars in a/ii* were taken into account. The 
best-fit values for each stellar property are obtained by 
computing the sum over all matches, weighted by their 
corresponding likelihood. Additionally, we account for 
the varying density of stars on each isochrone prescribed 
by the Initial Mass Function (IMF), by multiplying the 
weights by the number density of stars at each location 
as provided with the Y 2 isochrones. The IMF adopted 
is a power law with a Salpeter index. The effect of this 
latter weighting is generally small. 

The results for each system are listed in Table [3l where 
we give in addition to the mass and radius the theoretical 
values of log g^, luminosity L+, absolute visual magnitude 
My, and evolutionary age. In several cases our results 
differ slightly from those reported in recent discovery pa- 
pers that use the same atmospheric parameters adopted 
here and apply the a/i?* constraint essentially in the 
same way we have. The present study represents a slight 
improvement for those systems due to the application of 
weights, as described above. 

The relative errors of the stellar masses and radii de- 
termined here have median values of about 6% and 4%, 

4 In addition to the iron abundance, the enhancement of the a 
elements can also have a significant effect on the inferred stellar 
properties. About half of the transiting systems have at least one 
study in the literature reporting abundances for several of the a 
elements (Mg, Si, S, Ca, Ti, C, and O): HD 209458 (3 studies), 
TrES-1 (2 studies), WASP-f, XO-1, XO-2, and the five OGLE 
systems (one study each). In all cases the average enhancement 
is not significantly different from zero. We therefore assume here 
that it is zero for all systems. 
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Fig. 1. — Evolutionary trac ks for all host s tars except GJ 436 
from the Yonsei-Yale models of[Yi ct al. (2001]). The numbering of 
the systems follows that in the tables. The masses and metallicities 
adopted for the tracks, as well as the location of each star on the 
diagram, are from the best fit to the observations. 

respectively, but with wide ranges depending on the pre- 
cision of the observables (2% to 13% for <jM t /M+, and 
1.3% to 12% for a R jR+). While all of the stars in the 
current sample are hydrogen-burning stars, the degree of 
evolution within the main sequence varies considerably 
(see Figure[T]). Some systems such as TrES-3 (#11) and 
XO-1 (#15) are near the zero-age main sequence; others 
like HAT-P-4 (#20) have already lived for -90% of their 
main-sequence lifetime. The age is a critical ingredient 
for the theoretical modeling of the structure and evolu- 
tion of cxoplanets (see, e.g.. iBurrows et aTll2007| ). yet it 
is among the most difficult properties to determine for 
isolated main-sequence stars. For the host stars of tran- 
siting planets the evolutionary ages span the full range, 
as seen in Figure [2] The difficulty mentioned above is 
most evident for the less evolved objects (M + less than 
about 0.9 Mq), which are all seen to have large error 
bars that can reach the upper limit of 14 Gyr considered 
here. For those cases the nominal ages reported in this 
work should be used with caution. 

The surface gravities inferred from the models (log<?*; 
Table [3]) have formal uncertainties that are typically 
about 5 times smaller than those measured spectroscopi- 
cally (log g S p CC ; Table[T]). This reflects the strength of the 
constraint provided by a/R*. The values of log<7 spe c and 
log g+ are compared against each other in Figure [3J On 
average they agree quite well (the mean O—C difference 
is —0.027 dex), and the rms scatter of the differences 
is 0.15 dex although three of the systems present differ- 
ences larger than 0.2 dex. An illustration of the supe- 
rior constraint afforded by a/i?* is shown in Figure [4] for 
OGLE-TR-132 and WASP-2, two of the more dramatic 
examples. In neither case is the parallax known. 

3. COMPARISON WITH OTHER MODELS 

There is undoubtedly some systematic error introduced 
by imperfections in the stellar evolution models them- 
selves, but this type of error is difficult to evaluate. Ex- 
tensive comparisons between models and observations us- 



10 




1.0 1.2 
Stellar mass (M Q ) 

Fig. 2. — Evolutionary age versus mass for transiting planet host 
stars based on the models bv IYi et all 1(20011 ). GJ 436 is excluded 
(see text). The numbering of the systems is the same as in the 
tables. 




4.0 4.2 4.4 4.6 
Log g. from models 

Fig. 3. — Surface gravities f or the planet ho st stars inferred 
from stellar evolution models by lYi et al] 1120011 ) , compared with 
those measured spectroscopically. The dotted line represents the 
one-to-one relation, and the systems are numbered as in previous 
figures. 



ing double-lined eclipsing binaries with very accurately 
measured masses and radii have shown that the agree- 
ment with theory is in general very good, an d is within 
a few percent for solar-type stars (see, e.g.. lAndersenl 
[19911: IPols et al.l[T997t ILastennet fc Valls-GabaudH 2002) . 
As a simple test, we considered a second set of models by 
iGirardi et al.1 (|2000l ) that has often been used by other in- 
vestigators in the field of transiting planets. We used the 
isochrone interpolation tools provided on the web site at 
the Osservatorio Astronomico di Padova 5 to explore the 
agreement with the observed values of T e g, [Fe/H], and 

5 http: //stev . oapd. inaf . it/~lgirardi/cgi-bin/cmd . 



6 



Torres et al. 




8000 



6000 
T .„ (K) 



4000 8000 



6000 
T = „ (K) 



4000 



Fig. 4. — Measured properties of two extrasolar planet host stars 
displayed on plots analogous to the H-R diagram. The constraint 
on the location of the stars based on the surface gravities (log g S pec) 
and temperatures is shown in the top panels, and the same for a/R+ 
versus T c g is shown in the botto m panels. Isoch rones are from the 
series of evolutionary models by lYi et all l!2001fl for the measured 
metallicity in each case (Table [TJ , and are shown for ages of 1 to 
13 Gyr (left to right) in steps of 1 Gyr. The a/R^ values are seen 
to provide a much better handle on the stellar parameters (mass, 
radius, etc.). 

a/Ft*, and to infer the stellar mass and other properties 
in the same way as above. There are small differences 
in the physical assumptions between these models and 
those from the Y 2 series, but overall they are rather sim- 
ilar. For this application the Padova models are available 
for metallicities smaller than Z = 0.030 (corresponding 
to [Fe/H] = +0.20), which allows us to compare results 
for nine of the transiting systems. 6 Figure \5\ displays 
this comparison for the stellar masses and radii, showing 
that in all cases the results from both models agree to 
well within their uncertainties. 

Similar tes t s wer e carried out using the models of 
iBaraffe et al.l (|1998| ). These models enjoy widespread use 
for lower-mass stars and brown dwarfs, but they are also 
computed for masses as large as 1.4 Mq. They are avail- 
able for three different values of the mixing length param- 
eter qml over restricted ranges in mass and metallicity. 
The value that fits the observed properties of the Sun 
(against which all models are calibrated) is qml = 1.9. 
For cooler objects near the bottom of the main sequence 
the value used almost universally is «ml = 1.0. Here we 
adopt the former value, since our stars are all more mas- 
sive than about 0.8 M Q with the exception of GJ 436, 
which we discuss in more detail in §0 For the mass 
range of interest and for «ml = 1-9 the IBaraffe et al.l 
( 1998) models are publicly available only for solar metal- 
licity, so the comparison was limited to transiting sys- 
tems with [Fe/H] within about 0.1 dex of the Sun, and 
with [Fe/H] uncertainties that keep them within 0.2 dex 

6 Higher-metallicit y Padova models have been published by 
ISalasni ch et al. (2000), but those employ physical as sumptions dif- 
ferent enough that they cannot be merged with the Gira rdi et alj 
pOOOh models we use here. In addition, from a practical point of 
view, the use of these higher-metallicity models is not yet imple- 
mented in the online interpolation routine as of this writing. 




0.8 1.0 1.2 1.4 
M, from Y 2 models (M Q ) 




K 0.8 



0.8 1.0 1.2 1.4 
R, from Y 2 models (R ) 

Fig. 5. — Compariso n of the stella r masses and radii de- 
rived from the Y 2 m odels JVi et al. 2001) and the Padova models 
HGirard i et al. 200(J) for transiting planet hosts, showing the excel- 
lent agreement. The dotted lines represent the one-to-one relation, 
and the numbering of the systems is the same as in previous figures. 



of solar. Of the six systems in this range, one (OGLE- 
TR-113) gave no solution consistent with the observed 
values of T c g and a/R+. This is because the Baraffe 
models are limited to ages less th an about 10 Gy r for 
thi s mass range, where as both the lYi et al.l (|2001h and 
the lGirardi et al.l (|2000f ) models indicate an age of about 
13 Gyr for this star. Whether OGLE-TR-113 is truly 
this old is unknown, ft seems at least as likely that some 
of the other observational constraints are in error. The 
remaining 5 systems show excellent agreement with the 
results from both the Y 2 models and the Padova models, 
within the formal uncertainties. 

These tests of the stellar evolutionary models are obvi- 
ously not exhaustive, and it may well be the case that all 
of the sets of models we considered have some deficiencies 
in common. However, the general pattern of agreement 
does lend some degree of confidence to the results. We 
proceed under the assumption that the systematic errors 
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in these calculations are not the dominant source of error 
in the stellar parameters derived here (except perhaps for 
OGLE-TR-113, as noted above). 

4. ADDITIONAL OBSERVATIONAL CONSTRAINTS 

As a further test of the accuracy of our stellar radius 
determinations, in this section we consider the consis- 
tency check provided by the near infrared (NIR) surface 
brightness (SB) relations, which yield the angular diame- 
ter <f> of a star directly in terms of its apparent magnitude 
and color. Unlike the parallax, which would in principle 
yield the stellar luminosity but is known for only five of 
the brighter systems in the sample, <j> can be computed 
for all host stars from existing photometry. We use the 
empirical relation 

log 0sb = Cl (V-K)+c 2 -0.2K (4) 

derived bv lKervella et al.l (|2004[ ). in which the coefficients 
are a = 0.0755±0.0008 and c 2 = 0.5170±0.0017, the K- 
band magnitude is in the Johnson system, and </>sb is the 
limb-darkened value of the angular diameter expressed 
in milli-arc seconds. The relation represented by eq. ([4]) 
is extremely tight, with a scatter well under 1%. Near 
infrared magnitudes are available for all stars from the 
2MASS catalog and were t r ansfor med to the Johnson 
system following ICarpenterl (|200lD . The best available 
V magnitudes collected from the literature are listed in 
Table [U 

Angular diameters from our stellar evolution model- 
ing in § 12.31 can be derived for comparison with </>sb by 
making use of our theoretical radii and absolute visual 
magnitudes in Tabled along with the apparent V mag- 
nitudes, using 

mod = 9.3047i?*/10°- 2(v ^ Mv+5) . (5) 

With the stellar radius expressed in solar units, the nu- 
merical constant is such that </> mo d is in mas. Neither 
this equation nor the previous one take into account in- 
terstellar extinction, although eq. ([J| is actually quite in- 
sensitive to extinction. We discuss this below. 

The comparison between the angular diameters from 
eq. and eq. is shown in Table [5] for all systems 
except for GJ 436, for reasons to be described in the 
next section, and the OGLE stars, which are likely to 
be significantly affected by extinction since they lie sev- 
eral kpc away near the Galactic plane. The uncertainties 
listed include all contributions from the photometry and 
model-derived quantities, as well as the errors in the coef- 
ficients of eq. ||3J). The precision of 0sb is typically several 
times better than that of <fi mo d- A graphical comparison 
is shown in Figure with the one-to-one relation repre- 
sented with a diagonal line. The good agreement over the 
full range of an order of magnitude in <f> is an indication 
that our model radii from § 12.31 are accurate. 

One effect of extinction in this analysis would be to 
make the values based on eq. (JSJ) appear too small. Ex- 
amination of the differences shows a hint of this, partic- 
ularly among the fainter (more distant) stars: the av- 
erage difference ((0 mo d — 0sb)/^sb) is +0.5 ± 1.8% for 
the 6 objects brighter than V — 11, and —3.9 ± 1.6% 
for the others. Since reddening values for individual sys- 
tems are presently unknown, and cannot be determined 
accurately enough from the information at hand, we are 
unable to correct for this. We are equally unable to turn 
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Fig. 6. — Angular diameters 0sb computed from the n ear in- 
frared surface brightness relation of Kcrvclla ct al. (2004) com- 
pared against the values derived from our modeling, making use 
of the apparent V magnitudes for all host stars and ignoring ex- 
tinction. The dotted line represents the one-to-one relation, and 
the systems are numbered as in previous figures. 

the argument around and use the </>sb values as further 
constraints in our modeling. By trial and error we find 
that a uniform reddening value between E(B—V) = 0.02 
and 0.03 is sufficient to reduce the average discrepancy in 
the angular diameters for the fainter stars to zero. This 
modest amount of reddening is consistent with expecta- 
tions for objects that are between 100 and 500 pc away, 
as are these. 

In addition to the angular diameters, Table [4] lists the 
parallaxes 7r mo( j and corresponding distances we derive 
from V and My, ignoring extinction. Hipparcos par- 
allaxes 7Thip are given for comparison as well, for the 
objects that have them. 

5. THE CASE OF GJ 436 

As the only M dwarf among the currently known 
transiting planet host stars, GJ 436 presents a special 
challenge. The Y 2 stellar evolution models we used 
in all other cases are not intended for the lower main 
sequence, as they lack the proper non-grey model at- 
mosphere boundary conditions to the interior equations 
that have been shown to be criti cal for cool objects 
(see, e.g.. IChabrier fc Baraffill997| ). The Padova mod- 
els of lGirardi et"alT i|2000) have the same shortcoming, 
although both of these models are perfectly adequate 
for hotter stars. On t he other hand, the calculations 
bv iBaraffe et al.l ([1998) with ckml = 1.0 are specifically 
designed for low-mass stars, and have in fact been in- 
voked by other authors for this system as a rough con- 
sistency check. For the most part, the previous mass 
determinations for GJ 436 have relied upon empirical 
mass-luminosity relations, but even those relations have 
their problems. Previous radius estimates for GJ 436 
have often rested on the assumption of numerical equal- 
ity between M+ and R+ for M stars (jGillon et alj|2007at 
iDeming et afll2007t ). 

The importance of the GJ 436 system is undeniable, 
as it harbors the smallest transiting planet that has been 
found to date (comparable in size to Neptune). Fur- 
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thermore, it is the nearest transiting exoplanet, at a 
distance of only 10 pc. Given these facts, it is some- 
what surprising that until recently the m ass of the star 
(~0.4 Mfp) was only know n to about 10% (|Maness et al.1 
120071 iGillon eTa!ll2007H ). on a par with the worst of the 
determinations in Table [3] (that of WASP-2, with poorly 
determined spectroscopic parameters and more than 10 
times farther away) . Similar limitations hold for the stel- 
lar radius. As a result, our knowledge of the planetary 
parameters has suff ered. iTorresI (|2007l ) showed that the 
iBaraffe et al.l (|1998|) models in their original form do not 
yield satisfactory results for GJ 436 since different an- 
swers for M* and i?* are obtained depending on which 
constraints are used (including the absolute magnitude, 
since a reliable parallax measurement is available from 
Hipparcos). Moreover, some of the other predicted quan- 
tities are at odds with empirical determinations. These 
difficulties had not been pre viously e mpha sized. 

However, as described by ITorresI (|2007l ). reliable pa- 
rameters can still be obtained from the models by mod- 
ifying the theoretical radii and temperatures in such a 
way as to preserve the bolometric luminosities. This rec- 
ognizes the fact that i?* and T e g as predicted by theory 
have been shown to disagree with accurate measurements 
for M dwarfs in double-lined eclipsing binaries (see, e.g. , 
Pc^perljl997l; IClausen et all 119991 : iTorres fc Ribad 12002 



Ribad 120031 : Op cz- Morales k Ribas 2005 |), whereas the 



luminosities appear to be unaffected (Dclfoss e et al 
20001 ITorres fc Ribasl l2002t iRibad 12009: ITorres et al 



2006) . ITorresI ()2007l ) applied simultaneously the observa- 
tional constraints on T e g, a/i?*, the color index J—K, and 
the absolute if-band magnitude Mr-, and allowed the ra- 
dius/temperature adjustment factor to be a free parame- 
ter. In this way a self-consistent solution was achieved for 
all parameters, and in addition the radius/temperature 
factor showed excellent agreement with previous esti- 
mates for other M dwarfs. The stellar mass and radius 



are M* = 0.452 



+0.014 
^0.012 



Mq and R± = 0.464 



+0.009 „ 
-0.011 ft 



R n . We 



adopt these results here, thus placing GJ 436 on a similar 
footing as the other transiting systems in that the prop- 
erties of the host stars are all based on current stellar 
evolution models. 

The results for G J 436 are listed in Table [3] along with 
those of the other 22 stars. A few words of caution are 
in order. As pointed out by ITorresI |2007), the predicted 
absolute visual magnitude for this star is unreliable due 
to missing molecular opacity sources shortward of 1 /im 
in the model atmospheres that are used as boundary 
conditions in the stellar evolution calculations (see, e.g., 
IBaraffe et all 119981 : IPeTfosse et alJl2000T ). Also, because 
of the unevolved status of the star, the nominal age is 
probably meaningless since the observational constraints 
allow any age within the range of 1 -10 Gyr explored in 
the modeling effort of ITorresI <|2007l ) . Finally, this mod- 
eling effort assumed that the stellar metallicity is solar. 
For GJ 436 this is a valid assumption since the measured 
composition is [Fe/H] = —0.03 ± 0.20 (see Appendix), 
but the uncertainty in [Fe/H] was not accounted for due 
to the lack of proper models, and therefore the errors of 
the stellar properties in Table [3] may be underestimated. 

A careful r e ader m ay wonder whether our use of the 
IBaraffe et afl (|1998ft models in §[3] as a check on the 
results from the Y isochrones is contradicted by our 
remarks on the radius/temperature discrepancies men- 



tioned above. We do not believe so, because the mass 
regime of GJ 436 is very different. The host stars of the 
other transiting planets are typically more than twice as 
massive as GJ 436. 

6. PLANETARY PARAMETERS 

The combination of the stellar properties in Table [T] 
with the light curve parameters in Tabled along with the 
measured orbital periods and velocity semi- amplitudes, 
yields the homogeneous set of planet parameters pre- 
sented in Table [5] In addition to the planetary mass 
and radius, we have computed and tabulated the plan- 
etary surface gravity and mean density, as well as the 
orbital semimajor a xis and other usef u l char acteristics. 
As pointe d out by ISouthworth et al.l |2007). Winn et 
al. f2007). iBeattv et al.l l|2007| L and others, the surface 
gravity of a transiting planet can be derived without 
knowledge of the stellar mass or radius, using only the 
radial-velocity curve and the transit light curve. We take 
this opportunity t o correct the general expression for 
log<7 p presented bv lSozzetti et al.l (|2007h valid also for 
the case of eccentric orbits, which neglected to account 
for the projection factor implicit in the impact parameter 
b as derived from the light-curve fits: 



log 5p = -2.1383 - log P + logK* 



(6) 



5 log i- 



l-e 2 



a/R* 1 



•2 log 



' a/R* 
Rp/R± 



log(l 



The numerical constant is such that the gravity is in cgs 
units when P and K+ are expressed in their customary 
units of days and m s -1 . 

The first properties one generally wants to know about 
a transiting planet are its mass and radius. The large 
size measured for the first transiting planet discovered, 
HD 209458b, has been widely debated and still presents 
somewhat of a challenge to theory. It is now accom- 
panied by several other inflated planets, underlying our 
incomplete knowledge of the physics of these objects. An 
updated version of the now classical diagram of M p ver- 
sus R p is shown in Figure [3 in which five other examples 
are seen to be at least as large as HD 209458b (#3), or 
perhaps even larger, given the uncertainties. 

The fundamental physical properties of the known 
transiting planets cover a considerable range — more than 
two orders of magnitude in mass, a factor of nearly 5 
in radius, and a factor of 3 in orbital separation — and 
these are discussed in the next section. But the geomet- 
ric properties that determine the light-curve shapes are 
also varied. In Figure [5] we present a "portrait gallery" 
of all known transiting planets. Stars and planets are 
rendered to scale, emphasizing the wide range of stellar 
types probed by the photometric searches. The orbital 
geometries are indicated with solid horizontal lines rep- 
resenting an edge-on orientation, and the path of each 
planet across the stellar disk shown with dotted lines 
at the actual impact parameter. (Of course, the data do 
not distinguish between positive and negative impact pa- 
rameters.) The resulting light curves calculated for the 
V band are all shown with the same vertical (flux) and 
horizontal (time) scale, to facilitate comparison. Depths 
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Fig. 7. — Mass-radius diagram for all transiting planets (HAT-P- 
2 is off the scale, with M p = 8.7 Mj up ). Lines of constant density 
are shown. 

vary quite significantly (3 mmag in V for HD 149026b, 
26 mmag for OGLE-TR-113b; see Table E}, as do the 
overall shapes of the transit events, which in several cases 
are grazing enough to depart from the canonical profiles 
often depicted in the literature. 

7. DISCUSSION 

Many investigators have sought and claimed possible 
correlations between various stellar and planetary pa- 
rameters of transiting systems. In principle such correla- 
tions could lead to important insights into the formation, 
structure, and evolution of exoplanets. A primary mo- 
tivation for presenting a more complete, accurate, and 
homogeneous set of these parameters in this work was 
to facilitate such studies. The relatively large array of 
properties now available offers the opportunity to find 
new correlations, or to revisit old ones incorporating ad- 
ditional variables. While it is beyond the scope of the 
present work to investigate all possible correlations with 
statistical rigor, in this section we check on three of the 
most intriguing and potentially important relations that 
have been proposed. 

7.1. Planetary mass versus orbital period 

iMazeh et all (l2005h were the first to point out the ap- 
parent corre lation between M v and P for transiting plan- 
ets (see also lGaudi et al.H2005T) . The original suggestion 
was based on only 6 systems, but additional discoveries 
have generally supported the trend of decreasing mass 
with longer periods, although the scatter has also be- 
come larger. This is shown in Figure \§jp„. HAT-P-2b 
would be an extreme outlier in this plot; we have ex- 
cluded it because it is so much more massive than the 
other planets and may belong to a different category of 
planet (see § I7.2[I . Similarly, we have excluded GJ 436b 
because it is so much less massive than the others, and 
may be a rocky or rock- ice planet rather than a gas giant. 
A simple linear fit is shown for reference (dashed line). 

We also investigated the scatter in this relation, seek- 
ing any "third variable" that might correlate with the 



residuals. We seem to have found such a third vari- 
able: the metallicity of the host star. Figure^ displays 
the O — C residuals from the top panel as a function 
of [Fe/H], indicating a rather clear correlation (dashed 
line): AM p — (+0.152 ± 0.050) — (1.17 ± 0.23) x [Fe/H]. 
After removal of this trend, the relation between M p and 
period becomes tighter (Figure [SJ;). The scatter in the 
mass-period relation is reduced from 0.26 Mj up in the 
top panel to 0.17 Mj up in the bottom panel. It seems un- 
likely that this is a statistical fluke. However, the scatter 
is still larger than the formal observational uncertainties, 
suggesting that these three variables are not completely 
determinative. 

What might be the implications of this metallicity de- 
pendence? It has been proposed that the mass-period 
relation is related to the process by which close-in ex- 
oplanets migrated inward from their formation sites, or 
more specifically, to the mechanism that halts migration 
at orbital periods of a few days. The trend of larger 
masses at shorter orbital periods could suggest that the 
halting mechanism depends on mass, and larger planets 
are able to migrate further in. The dependence on metal- 
licity may then be interpreted to indicate that planets in 
metal-poor systems need to be more massive in order to 
migrate inward to the same orbital period as more metal- 
rich planets. In this context, the trend in Figure O) 
could be interpreted as evidence that the efficiency of 
the migration (or halting) mechanism is affected to some 
degree by the chemical composition. A dependence of 
migration on metallicity is in fact predicte d by some the- 
ories, and could arise, as pointed out by iSozzetti et al.l 
(2006), either from sl ower migration rates in metal-poor 
protoplanetary disks (|Livio fc Pringldl200l IBosd I2005h 
or through longer timescales for giant planet forma- 
tion around metal-poor stars, which would effectively 
reduce the efficiency of migration before th e disk dis- 
sipates dIdi~&Lin| [2001 [AlfoerFiFH3[200l). However, 
the above processes are more aimed at addressing the 
apparent lack of short-period planets amo ng very metal- 
poor stars claimed by some authors (see ISozzetti et al.l 
2006), whereas among the transiting planets it is not the 
lack of more metal-poor examples we are concerned with, 
but rather their different properties (such as mass) com- 
pared to metal-rich planets at the same orbital period. 
To give a quantitative example, we find that for a period 
of 2.5 days (near the average for known transiting sys- 
tems), the mass of a planet with [Fe/H] = —0.2 is ~40% 
larger than one with average metallicity ([Fe/H] = +0.13, 
excluding HAT-P-2b and GJ 436b), while the mass of a 
planet with [Fe/H] = +0.4 is about 30% smaller. This 
range of mctallicitics, from [Fc/H] = —0.2 to +0.4 (cover- 
ing a factor of 4 in metal enhancement), is approximately 
the full range observed. We note that the strength of the 
metallicity effect is modest rather than overwhelming, 
and this to o is in agreement wit h theoretical expecta- 
tions (e.g., EMoXPringli HH). The massive planet 
HAT-P-2b obviously does not conform to this trend of 
M p versus period, which may indicate some fundamen- 
tal difference either in its formation or migration. 

An alternative interpretation, also proposed by 
IMazeh e~a l. (2005), is that the M p versus P relation is 
more a reflection of survival requirements in close prox- 
imity to the star, due to thermal evaporation from the 
extreme UV flux. Close-in planets must be more mas- 
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Fig. 8. — "Portrait gallery" of transiting extrasolar planets. Star and planet sizes are shown to scale. The vertical and horizontal axes 
of all light curves are also on the same scale. Planet trajectories are shown with their measured impact parameters (dotted lines). Orbital 
periods are indicated in each panel. 
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Fig. 9. — (a) M p as a function of period for all transiting planets 
except HAT-P-2b (M p = 8.7 M Jup , P = 5.63 d) and GJ 436b 
(M p = 0.073 Mj up , P = 2.64 d); see text. The dashed line is a 
linear fit. (b) O — C residuals AM P from the linear fit in the top 
panel shown as a function of the metallicity of the host star. The 
dashed line is a linear fit. (c) Same as (a), with the dependence 
on [Fe/H] removed based on the fit in (b). The fitted line has the 
expression M p = (+1.70 ± 0.13) - (0.281 ± 0.044) X P. 
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Fig. 10. — Surface gravity versus orbital period for all transiting 
planets except HAT-P-2, which is off the scale (g p ~ 234 m s -2 ). 
A linear fit is shown. 

sivc to avoid ablation to the point of undetectability. 
The role of metallicity in this case w ould be through th e 
difference in the internal structure (|Santos et al J 12006). 
Metal-rich planets have been suggested to be more likely 
to develop rocky cores ([Pollack e t al. 1996; Guillo t~et al.l 
[200llBurrows et alj|2007l see also S ITS]) . If the presence 
of such a core somehow slows down o r prevents complete 
evaporation, as has been propos ed (|Baraffe et alj 12001 
iLecavelier des Etangs et al1l2004r ). survival at a given pe- 
riod would then have a dependence on metallicity. 

A diagram related to the one considered above is that 
of planetary surface gravity ve r sus o rbital period, first 
presented bv ISouthworth et"aT1 (|2007f ). Because g p does 
not require knowledge of the mass or radius of the star, it 
is in a sense a cleaner quantity that should be free from 
systematic errors in M* and i?* and is nearly independent 
of stellar evolution models. This not true of other bulk 
properties such as the mean planet density. An updated 
version of the g p versus P relation is shown in Figure flOl 
along with a linear fit. We examined the residuals from 
this linear fit for a possible correlation with stellar metal- 
licity, as we did earlier for the case of the M p versus P 
diagram, but we found none. Given that g p oc M p /R 2 , 
we note that the apparent influence of metallicity on the 
planetary radii must be contributing significantly to the 
scatter in the g p versus period relation. We discuss this 
further in §17.31 



7.2. Safronov number versus equilibrium temperature 

Recently, lHansen fc Barmanl (|2007l ) proposed a dis- 
tinction between two classes of hot Jupiters, based on 
a consideration of the Safronov number and the zero- 
albedo equilibrium temperature. The Safronov number 
is a measure of the abilit y of a planet to gravitation- 
ally scatter other bodies (Saf ronovl |1972| ). and is de- 
fined as 6 = ±{V csc /V ovh ) 2 = (a/R p )(Mp/M ir ), the ra- 
tio between the escape velocity and the orbital velocity 
squared. We assume that the zero-albedo equilibrium 
temperature scales as T eq = T e s(R*/2a) 1 / 2 (i.e., we as- 
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Fig. 11. — Diagram of Safronov number (0) versus equilibrium 
temperat ure (T eQ ) for t ransiting planets. Planet classes are labeled 
following Hansen & Barman (2007). A tentative dividing line at 
= 0.05 is indicated. 

sume that the heat redistribution factor / is common to 
all planets, in the absence of more complete knowledge). 

We l ist 9 and T ea for all tra nsiting planets in Ta- 
ble [3 lHansen fc Barmanl (|2007l ) pointed out a gap in 
the distribution of Safronov numbers, and defined Class I 
planets as those with 8 ~ 0.07 ± 0.01 and Class II as 
9 ~ 0.04 ± 0.01. They tentatively proposed also that 
these two categories have other distinguishing character- 
istics, such as a difference in the average temperature 
of the host stars, or the orbital separations . Upon the 
disco very of HAT-P-5b (jBakos et al.l l2007d) and HAT- 
P-6b (jNoves et al.ll2007l ). those authors pointed out that 
the distinction now seems less clear, as these two planets 
tend to fill the gap between Class I and Class II. 

The larger sample now available, and especially the 
more accurate and homogeneous set of properties pre- 
sented here, offers the opportunity to revisit the issue. 
Figure [Tl] shows an updated version of the 9-T oq dia- 
gram for transiting planets, which has some significant 
differences compar e d to t he original version. Following 
lHansen fe Barmanl (|2007f ) we have excluded the mas- 
sive planet HAT-P-2b, with a Safronov number so much 
larger than all the others (9 = 0.94) that it would seem 
to be in a different class altogether, as well as GJ 436b, 
which is of much lower mass and a presumably different 
composition. 

In our updated diagram the separation between Class I 
and Class II is still quite striking. The clustering of the 
Class II objects around the value 9 ~ 0.04 has tightened, 
if anything, and HAT-P-5b (#21) and HAT-P-6b (#22) 
do not encroach on the gap in Safronov numbers with 
Class I. (The Class II object with the lowest value of 9 
is OGLE-TR-ll lb f#6j.) Thus, the dicho tomy remains 
very suggestive. lHansen fc Barmanl (|2007l ) have argued 
that the principal distinction between the two classes is 
based on mass (planetary and/or stellar), and that plan- 
ets of Class II are, on average, less massive than those in 
Class I, and orbit stars that are typically more massive. 
While we agree with the latter part of this statement 
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Fig. 12. — Planetary mass as a function of equilibrium tem- 
perature for transiting planets. Class II planets (0 < 0.05) are 
represented with squares. A tentative dividing line is indicated. 

regarding the stellar masses, we find that the distribu- 
tion of planetary masses is indistinguishable between the 
two groups. We do confirm, however, that an updated 
diagram of M p versus T eq (see Figure [12]) shows planets 
in Class II to be systematically less massive for the same 
equilibri um temperature (level of i rradiation), as was also 
found by Han sen fc Barmanl (|2007f ) , so in this sense their 
general claim that the distinction has something to do 
with mass seems to be supported by the observations. 

An important characteristic that seems to be differ- 
ent in the two groups is the metallicity of the parent 
stars. This is illustrated in Figure [T31 Parent stars of 
Class II planets tend to be slightly more metal-rich. A 
Kolmogorov-Smirnov test of the two metallicity distribu- 
tions, which appear to be centered around [Fe/H] ~ 0.0 
for Class I and [Fe/H] ~ +0.2 for Class II, indicates 
only a 1.7% probability that they are drawn from the 
same parent population. There is perhaps a hint that 
the Safronov numbers for Class I planets show a decreas- 
ing trend with metallicity in Figure 1131 whereas the 9 
values for Class II planets are independent of [Fe/H]. 

7.3. Heavy element content versus stellar metallicity 

Transiting planet discoveries have challenged our the- 
oretical understanding of these objects from the very be- 
ginning. The inflated radius of HD 209458b, argued to 
be due to an overlooked internal heat source in the planet 
(see, e.g.. iFabrvckv et al"1l2007l and references therein), 
still poses a problem for modelers, and this first example 
is now joined by several other oversized planets indicat- 
ing it is not an exception (see §[6]). On the other end 
of the scale, HD 149026b was the first transiting giant 
planet found to have a radius si gnificantly smalle r than 
predicted by standard theories (jSato et "all 12005). The 
implication is that it must have a substantial fraction 
of heavy elements of perhaps ^70 M ffi (2/3 of its total 
mass), which is o ften assumed to be in a core. More re- 
cently HAT-P-3b ([Torres et al.ll2007t ) has also been found 
to be enhanced in heavy elements on the basis of its small 
size for the measured planet mass. In this case metals 



Improved parameters 



13 



0.08 



0.06 



0.04 - 



0.02 



-0.4 




120 



0.6 



Fig. 13. — Safronov number as a function of stellar metallicity 
for transiting planets. Class II planets (© < 0.05) are represented 
with squares. 

make up ~ 1/3 o f the t otal mass. 

iGuillot et all (|2006fl and also Burrows et ail (|2007l ) 
have found evidence that the heavy element content cor- 
relates with the metallicity of the parent star, a trend 
that was not anticipated by theory. These studies were 
based, respectively, on the 9 and 14 transiting planets 
known at the time. The larger sample now available 
warrants a second look at this possible correlation, for 
its potential importance for our understanding of planet 
formation. 

Here we have estimated the heavy eleme nt content Mz 
of eac h planet using the recent models bv iFortnev et al.l 
(|2007h . which include the effects of irradiation from the 
central star. The mass in heavy elements was calculated 
from the measured mass of the planet, its radius, the or- 
bital semimajor axis, and the age of the system as derived 
above from stellar evolution models. The uncertainty in 
these estimates is often very large, due mostly to errors 
in the measured values of R p and especially the age. In 
ten cases the models indicate no metal enhancement at 
all, and some of these are in fact "inflated" hot Jupiters. 
In a few other cas es only an upp e r limi t can be placed. 
For HAT-P-2b the IFortnev et at] (|2007l ) models as pub- 
lished do not provide a large enough range of core masses. 
The results for Mz are listed i n Tab le [5j A comparison 
with v alues from lGuillot et a l. (2006) and lBurrows et aTl 
(|2007h for the subset of planets in common indicates sig- 
nificant differences in some cases, either in Mz or its 
error. 

The mass in heavy elements is plotted against the stel- 
lar metallicity in Figure 1141 The overall trend is simi- 
lar to that pointed out by previous authors, in the sense 
that the upper envelope of the distribution appears to in- 
crease with [Fe/H]. It is natural to expect that a higher 
stellar metallicity implies a higher-metallicity protoplan- 
etary disk. In the context of the core-accretion model of 
planet formation, one would also naturally expect a more 
metal-rich disk to lead to more metal-rich planets, and 
in this sense the observed trend is in accordance with the 
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Fig. 14. — Heavy element content for transiting giant planets as 
a function of the metallicity of the host star. HAT-P-2 is excluded 
(see text). 

core- accretion theory. 

8. FINAL REMARKS 

Progress in understanding planet formation, structure, 
and evolution depends to a large extent on an accurate 
knowledge of their physical characteristics. This begins 
with an accurate knowledge of the properties of the par- 
ent stars. The main contributions of this work toward 
that goal are a re-analysis of all available transit light 
curves with the same methodology and the uniform ap- 
plication of the best possible observational constraints to 
infer the stellar mass and radius. This includes, in par- 
ticular, the constraint on the stellar density through the 
light-curve parameter a/R*. We expect the application 
of this technique to current and future ground-based or 
space-based transit searches to be highly beneficial, es- 
pecially when parallax information for the candidates is 
unavailable. 

Through the procedures described above we have ob- 
tained a more homogeneous set of stellar and planetary 
parameters than previously available, with error bars 
that are well understood and more appropriate when 
searching for patterns and correlations among the var- 
ious quantities. The results have enabled us to explore 
the role of metallicity in two of the more intriguing cor- 
relations between s tar and plane t prop erties: that of M p 
versus P, to which Ma zeh et alJ (|2005l ) first drew atten- 
tion, and tha t of the Safronov number versus equilibrium 
temperature (|Hansen fc Bar man 2007). We have also re- 
examined the correlation between the stellar metallicity 
and the heavy-element content of the planets. We find 
a clear influence of [Fe/H] in the M p versus P diagram, 
with the planets in more metal-poor systems being sys- 
tematically more massive at a given orbital period. This 
can be interpreted as evidence of a (mass and) metallic- 
ity dependence of the migration process. Alternatively 
it may be seen as indirect support for the correlation 
between core size and [Fe/H] along with the idea that 
the presence of such cores slows down or prevents com- 
plete evaporation of the planets in the extreme radiation 
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environments of these hot Jupiters. We find also that 
the improved parameters for the present sample of tran- 
siting planets support the recently proposed notion of 
at least two distinct clas ses of hot Jupiters in te rms of 
their Safronov numbers (|Hansen fc Barmanl 12007). Ad- 
ditionally, the average metallicity of Class II planets ap- 
pears to be ~0.2 dex higher than Class I. Thus, chemical 
composition adds another distinguishing characteristic to 
those proposed earlier. A further result of our work, 
which takes advantage of the larger sample of transiting 
planets now available, is the confirmation of the gen- 
eral trend of higher heavy-element content for planets 
with more m e tal-ric h ho st stars, originally ad vanced by 
iGuillot et all (|2006ft and lBurrows et all (|2007t ). 

We expect that the search for other significant corre- 
lations will be made easier by the better accuracy and 
homogeneous character of the properties of the present 
sample of transiting planets, and should lead to a deeper 
understanding of their nature. 
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APPENDIX 

SPECTROSCOPIC AND PHOTOMETRIC INFORMATION FROM THE LITERATURE 

Below we describe in detail our compilation of the information available in the literature for the atmospheric pa- 
rameters of all known transiting planets, as well as the sources for the light curves we have used, and for the velocity 
semi-amplitudes that determine the planet masses. For the atmospheric parameters (T e g, [Fe/H], log<?) we consider 
only spectroscopic studies, and disregard photometric determinations with few exceptions. The values adopted, which 
are listed in Table [TJ are the result of a careful combination of results for each planet rather than the selection of a 
particular favorite study, and we have therefore deemed it important to document our choices for future reference. 
For the spectroscopic parameters the goal of this effort has been to arrive at values that best represent the stellar 
properties based on current knowledge. Combined with our application of uniform procedures for deriving the stellar 
mass and radius based on state-of-the-art stellar evolution models, we hope that the overall accuracy of the planet 
parameters has been improved. 

1. HD 1490 26: Two co nsistent determinations of the effective temperature are available for this star. One is 
from iMasana et ail (|2006ft . who apply a semi-empirical method based on a spectral ener gy distribut i on fit (SEDF) 
calibrated against solar analogs, and the other is from a spectral synthesis analysis by Sato et aTl (12005ft carried 
out on high-resolution Subaru spectra using the SME package (Spectroscopy Made Easy. rValenti fc Piskunovlll996t 
iValenti fe Fischerl [2005ft . We adopt th e weig hted average of the two, with a generous uncertainty of 50 K. The iron 
abundance is adopted from ISato et ail J2005) , but with a more conservative e rror of 0.08 dex. Th e surface gravity is 
taken from the same source. The velocity semi-amplitude K is adopted from lButler et all (12006ft . We c onsider that 
study, based on 16 RV measurements with the Keck telescope, to supersede the one bv lSato et alj (|2005ft . which used 
fewer Keck velocities (6) combined with 4 velocities with the Subaru telescope, but obtained essentia lly the same result . 
Additional Keck velocity measurements (many obt ained during transit ) were reported recently by I Wolf et al.l (|2007ft . 
who also re-reduced the earlier Keck observations of lButler et al.l (|2006ft (for a total of 35) . Their study focused on the 
Rossiter-McLaughlin effect, and di d not present a va lue for K. The photometric data used here are the three Stomgren 
(b+y) /2 light curves presented by ISato et alj (|2005ft . t he Sloan q and r l ight curves presented by ICharbonneau et al.l 
(|2006ft . and the five (b+y)/2 light curves presented by I Winn et all (l2007eh. 

2. HD 189733: For T e g we adopt the weight ed average of the results bv lBouchv et al.l (|2005aft based on high-resolution 
CORALIE spectra, and lMasana et al.l (|2006ft . with an unc ertainty of 50 K Later papers by the first team quote values 
for T c ff (as well as log g and [Fe/H]) said to be taken from iBouchv et aT . (2005a), but which are not exactly the same 
for reasons that are unclear. We ado pt logg from the B ouchy work, as well as the metallicity, but with an increased 
[Fe/H] error of 0.08 dex. A study by iGrav et alj (|2003ft based on lower resolution spectra (1.8 A) gave a somewhat 
lower temperature (4939 K) along with a mu ch lower value for [Fe/ H] of —0.37; we do not consider those results here. 
The only published det ermination of K is bv lBouchv et alj (|2005af ). which we adopt. Additional RV measurements at 
Keck were obtained by Win n et al.l (|2006ft . mostly during transit for the purpose of studying the Rossiter-McLaughlin 
effect, but no value of K was reported. Ground-based light curves have been obtained in a variety of passbands 
since the discovery (Stromgren (b+y)/2, V, I, Sloan z). Hig h-quality observat ions with HST in a passband with an 
effective wavelength similar to the R band were reported by iPont et alj (|2007bft . These optical observations all suffer 
to some degree from the fact that the star is relatively active. Spott edness causes irregularities in the photom etry, 
which is very useful to establish the rotational period of the star (see IWinn et~all l2007ri iHenrv & Wind 12007ft . but 
interferes with the accurate determination of the light-curve parameters. By contrast. iKnutson et al.l (|2007ft obtained 
photometry with the Spitzer Space Telescope at the much longer wavelength of 8 /im, where spots have a negligible 
effect. Limb-darkening is virtually non-existent at this wavelength as well, removing another complication. While the 
ground-based, HST, and Spitzer light curve solutions all agree very well within their uncertainties, for the reasons 
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just described we have chosen in this case to rely only on the Spitzer data, which we have re-analyzed using the same 
methodolo gy applied to all the oth er transiting planets. A very similar (but less precise) radius ratio was reported 
recently bv lEhrenreich et all (|2007t ) based on additional Spitzer observations at 3.6 and 5.8 fim. 

3. HD 209458: Ten ind ependent determinations of T e g are available, which show excellent agreement: eight are 
from spectroscopic studies f Mazch ct al. 2000; Gonzalez et al. 2001; Mashon kina fc Ge hrcn 2001; Sadak ane et al.ll2002t 
iHeiter fc Lucid 2003: Sa ntos et al.ll200 4l : IValenti fc Fisch er 2005), and two others are based on photometry use the IRFM 
(|Rarmrez fc Melendezl I2004D and SEDF (|Masana et al J 120061) . We adopt the weighted average, with a conservative 
uncertainty of 50 K. Similarly, we adopt the weighted average of the nine [Fe/H] determinations from the spectroscopic 
sources, with an error of 0.05 dex. The log g value we use here is also the weighted average of the ei ght available 
measu reme nts. Four values for the semi-amplitude K have been published. The early determinations by Mazeh et al.l 
(2000) and I Henry et al.l (2000) were superseded by more recent studies based on many more RVs and refined analysis 
techniques by iNaef et al. ( 20041 ) (46 ELODIE and 141 CORALIE measurements over 5 years, not all published) and 
iButler et al.l (|2006[ ) (64 Keck measurements over 6 years) . We adopt the weig hted average of the results from the two 
new sources. Our light curve fits are based on the HST/STIS photometry of iBrown et alj (|2001l ). In this case, given 
the high signal-to-noise ratio of the data, we fitted for the quadratic limb-darkening coefficients along with all of the 
other parameters. 

4. OGLE-TR-10: A spectros copic study by l Santos et ail (|2006t ) used 4 co-added UVES spectra to derive T off , [Fe/H], 
and log g. An earlier study by IBouchv et al.l (|2005b[ ) used a subset of these spectra and somewhat different analysis 
techniques, and obtained a significantly hotter temperature and a higher metallicity. They also reported an unusually 
large surface gravity for the temperature they obtained (logg = 4.70, for T e g = 6220 K). Because these results may 
be affected by the st rong correlations often present between all these qua n tities, we have preferred not to use the 
Bouchv et alj (l2005bl) valu es here. Anoth e r rece nt study by iHolman et al.l (|2007aD supersedes previous analyses by 



Konacki et al.l ( 2003af ) and iKonacki et al.l (|2005h based on similar material, which, nevertheless gave essentially the 



same results with only slightly different errors. The Santos and Holman values for T D g (6075 ±86 K and 5800 ± 100 K, 
respectively) are somewhat discrepant (by 275 K, or 2.1ct). We note that the direction of the difference is consistent 
with the sign of the difference in the iron abundances of these two studies (Santos giving +0.28 ± 0.10, and Holman 
0.0 ± 0.2), given the typical correlation between these two quantities. Both teams have presented evidence supporting 
their T e g determinations (a good fit to the Ha line profile in the case of Santos et al., and agreement with the H/3 
profile as well as a consistency check using line-depth ratios, for Holman et al.). Under these circumstances we simply 
adopt intermediate values (arithmetic averages) for T g, [Fe/H], and log g, with increased uncer tainties of 130 K, 0.15 
dex, and 0.10 dex, respectively. Two consistent determinat ions of K have been r eported: one bv lBouchv et al.l ((2005b) 
based on 14 RVs from UVES/FLAMES, and the other by IKonacki et all (|2005). We adopt the latte r results, because 
they are based on a combination of 9 K eck measurements w ith the 14 RVs from lBouchv et al"1 (|2005bf ). The light-curve 
parameters are taken from the work of lPont et a l. (2007a), who combine the best photometry available. 

5. OGLE-TR-56: The adopted temperature, metalli city, and surface gravity for this star are the weighted average s 
of three fairly consistent spectroscopic determinations (|Konacki et al.ll2003bl : IBouchv et al.ll2005bl : ISantos et alJ l2006) . 
Conservative errors of 100 K, 0.10 dex, and 0.16 dex ar e assigned to t hose q uantities, respectively. Two determinations 
of the sem i-amplitude K are con sidered here: one by iTorres et al" (|2004[ ). who used 11 Keck RVs and whose study 
supersedes IKonacki et al.l (|2003bD . and one bv lBouchv et alj ( 2005c ). based on 8 UVES and 5 HARPS me asurements. 
We ado pt the weighted average of the two results. Our light curve fits are based on the V and R data of iPont et al.l 
(j2007al ). 

6. OGLE-TR-11 1: The atmospheric parameters for this system are adopted from the spect roscopic study by 
ISantos et al.l (|2006l ). based on 6 UVES spectra, which supersedes those in the discovery paper by IPont et al.l (|2004l ) 
that used only a subset of th e same spec t ra. T he errors seem realistic, so we adopt those as well. The only determi- 
n ation of K available is from lPont et alj (|2004l ). Our photometric solutions are based on the two /-band light curves 
of iWinn et all (p007l) . 

7. OGLE-TR-113 : Three spectroscopic determinations of T e g, [Fe/H], and logg give very consistent results 
(|Bouchv et al.l 12004 IKonacki et all 120041 : ISantos et all 120061 ). and we adopt the weighted average. The formal un- 
certainties are very small because of the (likely accidental) good agreement; to be conservative we have increased them 
to 75 K, 0.08 dex, and 0.22 dex, respectively. For K we have taken the weighted average from the first two sources 
above, which are based, respectiv ely, on 8 spectra from UVES and 7 from Keck. Our light curve fits use the R band 
photometry of iGillon et al.l (|2006f) . 

8. OGLE-TR-132: The temperature, metallicity, and surface gravity are t aken from IGillon et al.l (|2007al) . This 
high-resolution study, based on UVES spectra, supersedes an earlier one bv IBouchv et al.l (120041) . which had very 
large uncertainties. The velocity semi-amplitude is a dopted from the work of lMoutou et al.l (|2004l ). who obtained a K 
value some 20% larger than in the discovery paper of IBouchv et al.l (|2004l ) despite h aving used the same 5 original RV 
measurements. The difference is due solely to the improved ephemeris provided by iMoutou et al.l (|2004| ) based on a 
new, high-quality light curve gathered with the VLT, co mpared to th e origin al ephemeris from the OGLE photometry. 
Our light curve fits use the VLT i?-band photometry of IGillon et al.l (|2007aD . 

9. TrES-1: Four h i gh-resolution spectros c opic analyses of this s tar gi ving consisten t resul ts have been published by 
lAlonso et all (|2004t ). ISozzetti et all (|2004h . iLaughlin et all (l2005j ). and ISantos et all (|2006l ). We adopt the weighted 
averages of those studies. The metallicity from lAlonso et afl (|2004l ) was assigned here an uncertainty (not originally 
reported) of 0.2 dex based on previous experience with the same type of material. As before, we have chosen more 
conservative errors of 50 K for T g and 0.05 dex for [Fe/H], to account in the latter case for the fact that ISantos et all 
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(|2006f ) find a hint of small systematic differences with ISozzetti et all ( 2004'). The mean metallicity is close to solar: 
[Fe/H] = +0.02 ± 0.05. A claim has been made bv IStrassmeier k Ricd {2004) of a much lower metallicity for TrES-1 
of [Fe/H] = —0.6 (no error given). This is so discrepant compared to other four determin ations that it is most likely 
incorrect. For K, the only value appear ing in the literature is that of lAlonso et al.l (|2004l ). based on 8 Keck spectra, 
which we adopt. lLaughlin et al.l (|2005l ) made 5 additional RV measurements, also with Keck, and combined them 
with the Alonso velocities to improve the planet mass. However, no value for K was reported and it is not possible 
to recover i t accur ately from other published quantities. Our light curve fits use the Sloan z-band photometry of 
I Winn et alJ (|2007dD . 

10. TrES-2: A single source is available for the atmospheric parameters {Sozzetti et al.l [2007). and those values 
are adopted here. Although the published uncertainty in T e g is only 50 K, a number of checks on the temperature 
are presented that suggest t he value is accura t e, and thus we accept that error. The semi-amplitude K is adopted 
from the discovery pa per by lO'Donovan et al.l (|2006l ). Our photometric solutions rely on the z-band light curves of 
IHolman et all (|2007bf ). 

11. TrES-3: The atmospheric properties were taken from Sozzetti et al. ( in preparation), while th e semi-amplitude 
K and the light-curve parameters are adopted from the discovery paper bv lO'Donovan et al.l (|2007f ). 

12. TrES-4: The atmospheric properties were taken from Sozzetti et al. ( in preparation), while th e semi-amplitude 
K and the light-curve parameters are adopted fro m the discovery paper bv Mandus hev etafl (|2007l ). 

13. WASP-1: The spectroscopic investigation bvlStempels et al.l (|2007l ) based on NOT/FIES spectra supersedes the 
results in the discovery paper by ICollier Cameron et al.l (|2007t ).and in fact the new study claims that the SOPHIE 
spectra in the old one were compromised by scattered light and uncertain normalization. Despite this, the results 
for the atmospheric quantities are quite similar. We adopt the parameters from the new study, although with more 
conserv ative uncertainties of 75 K, 0. 08 dex, and 0.16 dex for T e g, [Fe/H], and log g. The only RV analysis available is 
that by ICollier Cameron et al.l (l2007f). based on 7 R V measurements with SOPHIE. Our light curve fits are based on 
the z-band photometry of ICharbonneau" eTaTl (p007h . 

14. WASP-2: For this sy stem there are no spectroscopic studies beyond the one in the discovery paper by 
ICollier Cameron et alJ (|2007l ). in which the uncertainties for T c g and log g are rather large. Nevertheless, the tem- 
perature appears to be accurate as shown by a consistency check obtained using the IRFM. The authors state that 
the abundances are not substantially different from solar, and they appear to adopt, according to the caption of their 
Table 3, a value of [Fe/H] = +0.1 ± 0.2. This may simply be a representative abundance, but it is a reasonable value 
since it is close to the average metallicity of the other transiting planets. We accept it for lack of a better determination, 
along with T c s and log g as reported. We take also the K value fr om this paper, which is bas ed on 9 RV measurements 
with SOPHIE. Our light curve fits u se the z-band photometry of ICharbonneau et al.l (|2007f ). 

15. XO-1: The discovery paper bv iMcCullough et al.l (|2006h presents the only spectroscopic study of this system. 
The atmospheric parameters are based on HJS spectra from the 2.7m telescope at McDonald Observatory, analyzed 
with SME. We adopt the temperature and metallicity from this source, with increased uncertainties of 75 K and 0.08 
dex. The log g value is adopted also as given there. The semi-amplitude K derived by these authors is based on 
10 RV measurement s from HET and the 2.7m (HJS). Our light curve fits are based on the z-band photometry of 
IHolman et all (|2006l ). 

16. XO-2: The values of T c ff, [Fe/H], and log g derived in the discovery paper bv lBurke et al.l (|2007l ) are adopted here 
as published, and are based on two spectra collected with the HET for this northern component of a 31" visual binary. 
The logarithmic iron abundance, +0.45 ± 0.02, is found to be very similar to that obtained for the visual companion, 
which is +0.47 + 0.02. We adopt conservative uncertainties of 80 K for T c g and 0.05 dex for [Fe/H]. The K value from 
this study is based on 10 RV measurements, and is adopted for our work. The light curve parameters adopted are also 
taken from the discovery paper. 

17 . HAT-P-1: A sing le determination of the atmospheric parameters has been published, in the discovery paper 
by iBakos et al.l {2007a) . It is based on an SME analysis o f a Keck spectrum. The star is the easterly and fainter 
component of an 11" visual binary known as ADS 16402. IBakos et al. (2007a) also analyzed the other component 
spectroscopically, and used both sets of parameters for this physically associated and coeval pair to better constrain 
the mass and radius using evolutionary models. In this process they noted some inconsistencies between the models 
and the observations (a steeper slope in the H-R diagram than predicted by theory), which advise caution in using 
the atmospheric parameters. We adopt the values as published, but with increased uncertainties of 120 K, 0.08 dex, 
and 0.1 5 dex for T c g, [Fe/H], and logg, respectively. The velocity semi-amplitude is also taken from B akos et all 
(2007a]), and i s based on 9 velociti es from Keck and 4 from Subaru. Our light curve solutions are based on the z-band 
photometry of IWinn et all (|2007ch . 

18. HAT-P-2: For convenience we use the HAT designation here, al though the s t ar is al so referred to as HD 147506. 
We adopt the atmospheric parameters from the discove ry paper bv IBakos et al.l {2007b), the only one presenting a 
detailed spectroscopic analysis. An independent study by lLoeillet et al.l ()2007| ) based on SOPHIE spectra did not yield 
better estimates, according to the authors, although it did con firm the [Fe/H ] value . The velocity semi-amplitude in the 
disco very paper has been superseded by the determination by lLoeillet et al.l {2007D . who included the original velocities 
from IBakos et al.l (|2007bf ) along with 63 new measurements obtained with SOPHIE. The latter were intended to study 
the Rossiter-McLaughlin effect, and were therefore taken mostly during transit. The authors modeled the complete 
velocity curve, and improved not only K b ut also the e c centric ity e and longitude of periastron lo. We adopt those 
elements in this work. In another study by IWinn et all (|2007a[) a total of 97 new spectra from Keck were obtained, 
also to investigate the Rossiter-McLaughlin phenomenon. The authors focussed mainly on that effect and did not 
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report a new value of K . The light curve parameters reported here are based on the z-band photometry of lBakos et al.l 
lj2007bl K 

19. HAT-P-3: Th e atmospheri c para meters are adopted from the only spectroscopic analysis so far, which is the 
discovery paper by iTorres et al.l (|2007f ) that is based on an SME analysis of a Keck spectrum. The uncertainties 
adopted for T e g, [Fe/H], and logg, which are 80 K, 0.08 dex, and 0.08 dex, are more conservative than the formal 
errors returned by SME, as described there. The value of K is taken from the same source, as are the light curve 
parameters. 

20. HAT-P-4: As in the previous ca se, the atmospheric parameters, the K value, and the light-curve parameters 
are adopted from the discovery paper (jKovacs et alJ 12007). The uncertainties in T e g and [Fe/H] are the same as for 
HAT-P-3. 

21. HAT-P-5: A ll atmospheric, spectroscopic, and light-curve parameters are taken from the discovery paper 
(|Bakos et al.ll2007d ). ' 

22. HAT-P-6: All parameters are taken from the discovery paper (|Noves et al.ll2007h . 

23. GJ 436: Onl y one valiant attem pt has been made to determine the atmospheric parameters of this M dwarf 
spectroscopically ()Maness et al.|[2007t ). It is based on both low- and high-resolution observations and a comparison 
with synthetic spectra. The difficulties of this kind of measurement are evidenced by the disagreements the authors 
obtain from their low- and high-resolution material, in both effective temperature (3500 K and 3200 K, respectively) 
and surface gravity (5.0 and 4.0 dex). They attribute this to shortcomings in our knowledge of the molecular opacities, 
mainly of TiO. They adopt a compromise value of T e ff = 3350 ± 300 K for a fixed surface gravity of log g = 5.0, and 
a metallicity set to the solar value. T he latter assumption is supported by a photometric determination of [Fe/H] 
= —0.03 ± 0.20 by iBonfils et al.l (|2005f) . based on the ab solute if -band magn itude and the V — K color. These are 
the values we adopt here. The velocity semi-amplitude by Ma ness et al.l (120071) uses 59 measurements from Keck, and 
supersedes the determination in the discovery paper by iButler et al.l (|2004h . T heir orbital model inc ludes a linear 
drift in the velocities presumably caused by a more distant orbiting companion. iDemorv et a l. (2007) use the same 
velocities, also solving for a residual linear trend, and they additionally incorporate an accurate timing measurement 
of the transit an d another of the sec ondary eclipse based on Spitzer observa tions. They r e port a K semi-amplitude 
similar to that of iManess et al.l (|2007l ). but give no uncertainty. We adopt the lManess et al.l (|2007h value here. For the 
light curve parameters we adopt a weighted ave r age of the values repo rted b y (or reconst r ucted f rom) the ground-based 
or Spitzer-based studies of iGillon et all (|2007bl ). iDeming et all (|2007f L and lGillon et alJ (|2007d ): a/R* = 13.34± 0.58, 
R p /R* = 0.0834 ± 0.0007, and b = 0.848 ± 0.010. 
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TABLE 1 

Adopted atmospheric properties and orbital semi-amplitudes of the host stars. 
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Note. — Sec Appendix for the sources of these determinations. Apparent magnitudes for the OGLE stars are in 
the / band rather than V , and have estimated uncertainties of 0.05 mag. 

a These spcctroscopically determined surface gravities arc generally less reliable than those derived from our stellar 
evolution modeling (see Table [3j, and arc included here only for completeness.* 5 The large photometric uncertainty 
in V for such a bright star is due to the variability of the object (sec, e.g., Winn ct al. 2007b). c The orbit of the 
planet is not circular, and has e — 0.5170^q'qq^q and uj — 189. I^q'^ degrees dLo cillct ct al. 2007).^ The orbit of 
the planet is not circular, and has e — 0.14 ± 0.01 and u) — 350 degrees (Dcmory ct alj 2 007T . No uncertainty was 
reported for uj. The orbital fit includes a linear trend presumably due to an outer orbiting companion. 
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TABLE 2 

Light curve parameters for transiting planet systems. 







Period 








i 




Depth in V 


# 


Name 


(days) 


R P /R* 


a/R* 


b = a cos i/ R+ 


(dog) 


(g cm 3 ) 


(mmag) 



1 


HD 149026 


2 


.87598 


0491 +°- 0018 
u.u £ tyi_Q poos 


7 ii+0.03 
' ■ 1± -0.81 


0.00±g;j$ 


90.0±g;g 


0.8221°;°" 


3.0 


2 


HD 189733 


2 


.218573 


0.154631™ 


o 01 +0-06 
°- 8i -0.06 


0.680t°;g» 


85.58±g;g| 


2-629lg° 5 5 4 


25.3 


3 


HD 209458 


3 


.524746 


0.12086l° ™ 




o-eo7±g : g« 


86.7llg 95 5 


1 024+ 014 
u — 0.014 


17.0 


■1 


OGLE-TR-10 


3 


.101278 


H0+ - 002 
u.liu_ 002 


8-07±g;** 


0-Ool°; 5 00 


90.ol°;° 


1 03+ - 18 
±.uo_ 24 


15.2 


5 


OGLE-TR-56 


1 


.2119189 


1027+ 0019 
U.±UZ(_o 0019 


3.74l° ;i| 


o.8l7tg;S! 


77.60±°; 91 


674+ 11 
u '°'*— 0.078 


10.0 


6 


OGLE-TR-111 . . . 


4 


.01610 


129Q+ 0010 
u.±zyy_Q 0013 


12.09i°;« 


0.360t°;° 7 8 


QQ occ + 0.83 
o».^0_o 3 


2 07+ 24 
z-u ' —0.22 


20.4 


7 


OGLE-TR-113 . . . 


1 


.4324752 


1450+ 0016 

U.110U_q .0005 


6.38±°:S? 


0.24l5;3 


87.80lJ'g 2 


9 qqc+0.034 
i,J3u _ 0.29 


26.0 


8 


OGLE-TR-132 . . . 


1 


.689868 


0932+ 0015 
u.uyo^_Q 0020 




0.560lg : ° 46 


83.41 2 ; 9 


72+ - 29 
U- ' —0.13 


10.0 


9 


TrES-1 


3 


.030065 


fl 1QC7S+0. 00035 
U.IOO i 0_q .00032 


10.52±°;°| 


00+ 19 
U.UU_o oo 


90.ol°;° 


2.400l°;°j 4 


23.2 


10 


TrES-2 


2 


.47063 


n iokq+0.0010 

u.izoo_ ooiO 


7.62±°;S 


U.8540_ 0062 


83.57I°; 14 


i o 70 +0.061 
l.J(Z_o 059 


13.6 


11 


TrES-3 


1 


.30619 


1660+ 0024 
U.lODU_ 0024 


fi nfi+ 010 

°- UD -0.10 


n OO77+0-0097 
/_ .oo97 


82.15I°; 21 


2.47I°; 12 


24.6 


12 


TrES-4 


3 


.553945 


09903+ 00088 

u.umuo_ o.()()()88 




7^+0.015 
U- '00_q .015 


°2.»l_o 33 


0.328lg° 22 


9.8 


13 


WASP-1 


2 


.519961 


1025+ 0007 
u.iuzo_ 0007 


5.691°;°? 


00+ 27 
U.UU_ go 


90.ot°;° 


n =40+0.009 
u.o^y_ 59 


13.2 


14 


WASP-2 


2 


.152226 


1S10+ 0013 
U.1J1U_ 0013 


7.95«-j 2 


n 724+0-017 
U- '^*_0.028 


84.8ll°-J 7 


2 05+°- 26 


17.7 


15 


XO-1 


3 


.941534 


1326+°°°° 4 


11.55lg;S 


0.240l°;°f 


oo O1+0.70 
S8.51_o 30 


1.877l° : °i 5 


21.8 


16 


XO-2 


2 


.615857 


insq5+ 00090 

U.1UJ90_ 00085 


« 9 + - 1 
8 -^-0.2 


0.158lo.085 


88.90l°; 7 ° 


i t; 9C :+0.057 


13.5 


IT 


HAT-P-1 


4 


.46543 


1124+ - 0007 
u - iiz4 -0.0007 


1 47+0-23 
iu - 4 ' -0.23 


71 2+ - 017 
u -' iz -0.017 


86.llt8-.il 


1 nRQ+ 074 


13.1 


18 


HAT-P-2 


5 


.633320 


0.0b840_ 00073 


10 12+ ' 03 


00+ 44 
u - uu -o.oo 


90.01°;° 


0.618lHe 6 


5.9 


19 


HAT-P-3 


2 


.899703 


1109+ 0025 

u.liu»_ 0022 


iu.oy_ g4 


51+ - 11 
U - Oi -0.13 


87.24lg;« 9 


2 67+ - 54 
z - D '-0.59 


14.3 


20 


HAT-P-4 


3 


.056536 


n ns9nn+ 00044 
0.08200_ 00044 


6 04+ 03 


o oi +0 - 23 

u - u± -0.01 


89.9ll°;° 9 


446+ 007 


8.5 


21 


HAT-P-5 


2 


.788491 


1106+ 0006 

U.11UD_ 0006 


'■ OU -0.19 


425+ 045 


86.75lg' 44 


, n97 +0.081 
1 - uz ' -0.076 


14.6 


22 


HAT-P-6 


3 


.852985 


n ncm8+ 00053 

u.uaooo_ .00053 


7 6Q + °' 22 
' - OM -0.22 


602+ 030 
U.DUZ_ 030 


85-5ilg;f 5 


U.b8U_ 04g 


9.7 


23 


GJ 436 


2 


.64385 


0.0834_ Q 0007 


1 3 -34+0-58 


n S4S+ 010 
0.848_o.oio 


86.36lgi° 


6 43+°- 87 
0.80 


6.5 



a Stellar mean density computed directly from a/R* and P using p k — -^(a/R*) 3 (see lWinn eTaUl2007el) . ignoring the very small 
corrective term due to the planet. 
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TABLE 3 

Inferred properties of the host stars from stellar evolution models. 



# 


Name 






log 
(cgs) 


£* 


(mag) 


Age 
(Gyr) 


1 


HD 149026 


1 294+?,'*" 


1.368_Q.og 3 


4 278+°'° 4 t 
' ° — 0.063 


° — 0.35 


3 80+°i® 

— 0.24 


1 9+°i 

J — 0.9 


2 


HD 189733 


806+°° 4 f 

u-ovvu — 0.048 


756+° °}! 

u * ' ,JXJ — 0.018 


4 587+°'°" 

uu ' —0.015 


n oq-, +0.029 
u,00i - 0.028 


6 20+° " 

U *^ U — 0.10 


6. 8+44 


3 


HD 209458 


i ii9+° °22 


1 155+°°" 

i.ivJJ_0 016 


4 361+°'°°! 


1 622+?,' 09 , 7 


4 273+°'°J§ 

-t.^ < o_ 0.069 


3 1+ ( ''S 

t, - i — 0.7 


4 


OGLE-TR-10 


1 14+ - 1 !! 


1 17+°'" 

-0.11 


4 358+°'°° 4 

4,JJO -0 .082 


1 54+2'5? 
Xi J -0.38 


4 34+°'?? 
-0.35 


3 2+ 4 '° 


5 


OGLE-TR-56 


0.078 


1 363+SSI? 

1 - JIW - 0.086 


4 258+?,' , 43 , 

^•^°° — 0.043 


2 24+°' 4 ? 

0.40 


3 92+S'?? 
^•^ — 0.23 


3 2+;° 

J -^— 1.3 


6 


OGLE-TR-111 . . . 


852+°'°g 

0.052 


831+°'°^ 

u * oa — 0.040 


4 529+°'° 3 * 

t ±-<->^_0.042 


401+8'SI? 

w w — O.Obl 


5 99+° ?° 


8 8+5'? 

° ° — 6.6 


7 


OGLE-TR-113 . . . 


779+°'°" 

o. i i ^_o.015 


774+°'°?° 

' 1 —0.011 


4 552+°'°°? 


296+° °?i 

w *^ u — 0.018 


6 400+°'°°° 


13 2+°'5 

— 2 . 4 


8 


OGLE-TR-132 . . . 


1 305+Pv°S 

iilJUvJ — 0.067 


1 32+°" 


4 313+°'°^ 

^• J - LJ -0 .090 


2 35+°'I° 

JJ -0.49 


3 84+°'?° 

° ° — 0.33 


i.i 


9 


TrES-1 


878+°°? 8 

°-0.040 


807+° °" 

u-oui _Q.oi6 


4 567+° °?? 

t-oui -0.015 


438+°'° 3 ,5 

u,4JO - 0.033 


5 847+°'°°! 

o.o'ti _Q.095 


3 7+ 3 ' 4 

J - ' -2.8 


10 


TrES-2 


983+S-^? 

u -" OJ - 0.063 


1 003+ , . 33 , 

1,ulJJ -0 .033 


4 427+°'°!? 
-0.021 


1 06+°' 10 


4 77 +o.i2 

1 ' -0.11 


5 0+?' 7 

J - u -2.1 


11 


TrES-3 


915+°-°^ 


812+° °" 

uo — 0.025 


4 581+° °!! 

J " -0.012 


592+2 °^ 


5 46+°'" 

^•^U— 0.13 


6+?? 


12 


TrES-4 


1 394+22S? 

i.o»*i_0 056 


1 816+°'*5 

ljOlu -0 .062 


4 064+°'°?! 


4 39+°,'^ 


3 i6+°' 13 

o.±u_ () 14 


2 9 +t '' 4 


13 


WASP-1 


1 301+°'° 4 ° 


1 517+S'SS 

i -0.045 


4 190+°°?° 

^• X£ " J _0.022 


2 88+ ,' 3 , 6 , 
^■ oo -0.30 


3 63+° 13 


3 0+°'° 


14 


WASP-2 


89+ ' 1 ? 


840+°'°°? 


4 537+°'°^ 

■i.uji -0.046 


47+° }° 

u '^' -0.13 


5 78+°' 4 ° 

u ' ' °— 0.36 


5 6+!' 4 

°- u -5.6 


15 


XO-1 


1 027+°°" 


934+° °" 


4 509+°,'°,l? 


86+2'}? 

u - ou — 0.10 


5 02+°' I 4 

°* u — 0.16 


1 0+^'i 


16 


XO-2 


974+° °^ 

ual -0.034 


971+°'°?I 

u a ' —0.026 


4 452+°'°?° 
°— 0.022 


689+2'°?? 

u,llra -0.074 


5 33+°' ; 4 


5 8+?'l 

° ° — 2.3 


17 


HAT-P-1 


. ,,0+0.075 
i.ioo_ 07g 


i ioc+0.048 
i.ijo_ 04g 


4 382+ ' 027 
^.joz_ 030 


1 48+ - 25 
1 ' 4S -0.25 


4 S8+ ' 21 
4 - J5 -0.19 


n 7 +2.5 
^•'-2.0 


18 


HAT-P-2 


1 Qns+°' 088 

1 - ouo -0 .078 


1 506+ 013 


4 iqq+0-043 
4.iyy_ 053 


o ori+0.82 
°' zu -0.58 


3 50+°' 24 
o.ou_ 27 


2 6+°' 8 


19 


HAT-P-3 


n Q9S+ 044 
0.928_ 0S4 


n oqq+0.034 
u.ooo_ 044 


4 564+ 032 
^•ODt_0.032 


u.44y_ 064 


5 84+ ' 21 
°.° 4 -0.19 


1-5151 


20 


HAT-P-4 


1.248l°;°2° 


1 59fi+°° 6 ° 
i.oyo_ 075 


- 197 +0.019 
4 - lz '-0.027 


2 ?n +0.37 
z - ,u -0.36 


o 74+0.17 
^■'^-o.ie 


4-6l 2 ;8 


21 


HAT-P-5 


. .,. 7 +0.043 
± - lo '-0.081 


1 165+ ' 046 
i.ioo_ 052 


4.368t°'° 2 f 


1 54+0-24 
± '° yi -0.22 


4 OO+0.19 


2 6+ 2 ' 1 


22 


HAT-P-6 


1 290+ 064 


I 4KQ+0.069 

±.tUO_Q Q 63 


- 91s +0.027 
i.^io_ 030 


3-5918:18 


3 36+°' 16 
°- ou -0.16 


9 O+0.5 
J -0.6 


23 


GJ 436 a 


4tS2+°- 014 


4fi4+° °° 9 


4.84318:81? 


0.0260l8°8 1 i 4 7 


in 244+ 082 

1U.Z44_q 0g2 


6 0+ 4 '° 
°- u -5.0 



a The value of My is biased due to missing opacity sources in the models (see text), and the nominal age is only 
indicative since it is essentially unconstrained for this unevolved M dwarf. 
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TABLE 4 

Angular diameters and parallaxes for transiting planet hosts. 







0mod 








Distance 


7I - mod 


^HIP 


# 


Name 


(mas) 




(mas) 




(PC) 


(mas) 


(mas) 


1 


HD 149026 


0.171 ± 0.021 





1747 ± 


0021 


74.4 ± 7.2 


13.43 ± 1.29 


12.68 ± 0.70 


2 


HD 189733 


0.357 ± 0.020 





3609 ± 


0055 


19.7 ± 1.0 


50.70 ± 2.55 


51.94 ± 0.87 


3 


HD 209458 


0.2269 ± 0.0083 





2211 ± 


0038 


47.4 ± 1.6 


21.12 ± 0.72 


21.24 ± 1.00 


1 


OGLE-TR-10 
















5 


OGLE-TR-56 
















6 


OGLE-TR-111 
















7 

8 


OGLE-TR-113 
OGLE-TR-132 
















9 


TrES-1 


0.0493 ± 0.0024 





0487 ± 


0006 


152.3 ± 6.7 


6.57 ± 0.29 




10 


TrES-2 


0.0439 ± 0.0028 





0451 ± 


0007 


213 ± 11 


4.70 ± 0.25 




11 


TrES-3 


0.0309 ± 0.0018 





0339 ± 


0004 


245 ± 13 


4.08 ± 0.22 




12 


TrES-4 


0.0349 ± 0.0025 





0352 ± 


0005 


485 ± 31 


2.06 ± 0.13 




13 


WASP-1 


0.0346 ± 0.0026 





0359 ± 


0005 


408 ± 27 


2.45 ± 0.16 




14 


WASP-2 


0.0471 ± 0.0096 





0560 ± 


0013 


166 ± 30 


6.03 ± 1.10 




15 


XO-1 


0.0506 ± 0.0041 





0531 ± 


0007 


172 ± 12 


5.82 ± 0.42 




16 


XO-2 


0.0610 ± 0.0045 





0609 ± 


0009 


148 ± 10 


6.75 ± 0.46 




17 


HAT-P-1 


0.0680 ± 0.0071 





0699 ± 


0011 


155 ± 15 


6.44 ± 0.61 




18 


HAT-P-2 


0.127 ± 0.018 





1169 ± 


0016 


110 ± 13 


9.07 ± 1.06 


7.39 ± 0.88 


19 


HAT-P-3 


0.0559 ± 0.0059 





0594 ± 


0011 


139 ± 13 


7.21 ± 0.68 




20 


HAT-P-4 


0.0474 ± 0.0042 





0457 ± 


0007 


314 ± 24 


3.19 ± 0.24 




21 


HAT-P-5 


0.0317 ± 0.0031 





0334 ± 


0005 


342 ± 30 


2.93 ± 0.26 




22 


HAT-P-6 


0.0522 ± 0.0046 





0534 ± 


0008 


261 ± 20 


3.83 ± 0.29 




23 


GJ 436 














97.73 ± 2.27 


Note. — </> m od and 


i>SB are the angular diameters derived using eq. |5j and eq 


till . respectively. 


The distance and 



parallax are computed from My (Table [3j and V (Table ignoring extinction. 



TABLE 5 

Properties of extrasolar transiting planets. 



M p R p logg p p p a T cq M z 

_# Name (M Jup ) (R Jup ) (cgs) (g cmj (AU) 9 (K) (M ffi ) 

1 HD 149026 0.3591°°^ 0.654«;^ 3.360t° °|| 1-59±°;H 0.04313±°;°°°i 0.0384t°;°°| 8 7 1634^ 68^° 

2 HD 189733 1.144j:°;°^ 1.1381°;°^ 3.340l°°^ 0.9631°;°^ 0.03099to.oooe3 °- 0772 -om)27 120lt^ 0-35 

3 HD 209458 0.685^ °^ 1.359±g;°i| 2.963±°;«°| 0.338±g;°12 0.047071°;°°°*? 0.04234t°;^8 1449 +12 „ 

4 OGLE-TR-10 0.62±°;J* 1.25±g:^ 3.00tg;^ 0.40^'^ 0.0434^°;°°^ 0.0386±g;gg|§ 1481+^ 0-25 

5 OGLE-TR-56 l-39±g;l? 1.363±°;°gg 3.2641°;°* 0.68± {l 0.02383l°;°oosf 0.0393±°;°°jg 22 1 2 ±63 

6 OGLE-TR-111 . . . 0.55t°;J° 1.051±°;°g 3.088l°;°92 0.59±°;lf 0.04689«;°° )( 1 J ° 7 0.0571°;™ 1025^ W±™ 

7 OGLE-TR-113 . . . 1.26±g;ig 1.093t°;°g 3.419±°;°g l-^to.le 0.022891°;°°°^ 0.0677l°;°os? 25 ±25 

8 OGLE-TR-132 . . . 1.18±°;£ 1.20^ 3.276l°;J° 4 0.85l°;g 0.03035l°;°0053 0.0440±°;°°°| 20131^ 5tf 

9 TrES-1 0.752±°;°£ l^l ,^ 3.220t°;°£ 0.769t°;°£ 0.03925±g : ggg« 0.06341°;°°* 1140+3 20^° 

10 TrES-2 l-20 ™ l-^-Hfi 3.298±g;gig 0.813±g:gH 0.03558±g : ggg™ 0.0709+°;°™ 1498±1? 

11 TrES-3 1.938±g;gg 1.312+°;°g 3.452+°;™ 1-065+°;^ 0.02272±g;gggJ5 0.0738+°;°™ 1623 +26 _ 30 

12 TrES-4 0.920+°;° 7 3 2 1.75ll°;gg 2.872+°;°^ 0.213+°;°f° 0.05092+°;°°°™ 0.0384+°;°°g° 178511 

13 WASP-1 0.918l°;° 9 90 1-514+°;°!? 3.010+°;° 44 . 0.328+°;°« 0.03957+°;°™g 0.0374+°;°°3 7 isil+J* 

14 WASP-2 0.915±g;ggg l-071±g;g|° 3.287+°;°f 3 0.92+°;£ 0.03138±g;ggit2 0.0596±g;gg« 1304±g* 20+?* 

15 XO-1 0.918t°;° 8 78 l-206+°;°2 3.211±S;°£ 0.650+°;™ 0.04928+"° 0.0744+°;°°« 1196+? 3 0-20 

16 XO-2 0.566±g;g» 0.983+°;™ 3.168±g;g« 0.741+°;^ 0.036841°;°^ 0.0438±g;ggg 1319+** 30+» 

17 HAT-P-1 0.532+°;° 3 : ° l-242±g;gg 2.931±g;g« 0.345+°;™ 0.0553±g;gg» 0.04181°;°°?° 1306+ 3 3 ° 

18 HAT-P-2 8.72+°;?,° 1.003±g;g|| 4.370t°;° 7 l 0.0679±g;ggiS 0.941+°;g?| 1398±|i 

19 HAT-P-3 0.596±°;°26 0.899i°;°« 3.310+°;g 7 ?, l-02±g;i2 0.03882±°;°°° 7 ° 0.0585+°;°° 4 4 1127+. 49 , 70±fg 

20 HAT-P-4 0.671+°;g 44 1.274+g;° 4 g 3.020±g;gJ| 0.403l°;°* 0.04438t°;°o?5 1 0.03781°;°°^ 1686+2,° 

21 HAT-P-5 l-06±g;ii 1.254+°;™ 3.219±g;g|I 0.66+°;^ 0.04071±g;ggg« 0.0591+°;ggg 4 153911 

22 HAT-P-6 l-059±g : ggi 1.330+°;°g 3.171+°;°f 0.559+°;° 8 7 g 0.05237+°;°°™ 0.0646+g;°° 3 3 ? 1675+1? 

23 GJ 436 0.0729+°;°™ 0.3767+°°™ 3.107+°;g 4 ° 1.69+g;« 0.02872+°;°°™ 0.02461°;°°" 649+°", 

Note. — O is the Safronov number, T eq the zero-albedo equilibrium temperature ignoring the energy redistribution factor, and Mz is the planet heavy 
clement content. 
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